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ABSTRACT 



We propose an indirect approach to the generation of a two-dimensional quasiperi- 
odic (QP) pattern in convection and similar nonlinear dissipative systems where a direct 
generation of stable uniform QP planforms is not possible. An eightfold QP pattern can 
be created as a broad transient layer between two domains filled by square cells (SC) 
oriented under the angle of 45 degrees relative to each other. A simplest particular type 
of the transient layer is considered in detail. The structure of the pattern is described 
in terms of a system of coupled real Ginzburg-Landau (GL) equations, which are solved 
by means of combined numerical and analytical methods. It is found that the transient 
"quasicrystallic" pattern exists exactly in a parametric region in which the uniform SC 
pattern is stable. In fact, the transient layer consists of two different sublayers, with a 
narrow additional one between them. The width of one sublayer (which locally looks like 
the eightfold QP pattern) is large, while the other sublayer (that seems like a pattern 
having a quasiperiodicity only in one spatial direction) has a width ~ 1. Similarly, a 
broad stripe of a twelvefold QP pattern can be generated as a transient region between 
two domains of hexagonal cells oriented at the angle of 30 degrees. 
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1 Introduction 



Stable quasiperiodic (QP) planforms have been first discovered as equilibrium patterns 
in metallic alloys M. The simplest and most important types of the QP planforms are 
the ten-, eight-, and twelvefold ones (note an argument according to which only these 
types of the two-dimensional quasicrystals may occur in physical systems [0]). A 2n-fold 
pattern is a superposition of n > 4 spatial harmonics based on a set of equal-length wave 
vectors kj with the equal angles n/n between kj and k J+1 (more general aperiodic patterns 
with unequal wave vector lengths and/or unequal angles between the vectors may also 
exist, but they have probably never been considered, except for the case of unequilateral 
hexagons with n = 3 j3|). 

Then, it was predicted that the 2n-fold structures of different particular types may 
exist as dynamical planforms of a stationary form in nonequilibrium systems, such as 
thermal convection and the like [|J. Later, the predicted pattern (a twelvefold one, in 
particular) was indeed experimentally observed, first, as a stable dynamical nonequi- 
librium planform in the Faraday ripples || (a large-aspect-ratio liquid layer subject to 
high-frequency shaking in the vertical direction), and recently in an optical cell filled with 
Na vapors ||. Nevertheless, experimental generation of dynamical QP patterns is far 
from being straightforward; in particular, in the work [g] it was necessary to shake the 
liquid layer with a two-frequency quasiperiodic force at a specially selected ratio between 
the two frequencies. These experimental difficulties are related to the general fact that, 
in terms of the corresponding coupled Ginzburg-Landau (GL) equations for amplitudes 
of the spatial modes, a superposition of which gives rise to the pattern, a formal solu- 
tion for the QP planforms always exists, but it may be stable only in a relatively narrow 
parametric region |4j]. 

An objective of this work is to put forward a much easier possibility of generating eight- 
and twelvefold QP planforms of a finite but large size (a stripe) in "normal" systems, where 
direct generation of the quasicrystallic patterns does not seem feasible. To this end, we 
notice that the set of the wave vectors on which the eight- or twelvefold QP is based may 
be regarded, in an obvious way, as a superposition of two half-sets of 2 or 3 vectors, see 
Fig. [T] below. Each half-set, in turn, may give rise to a usual periodic pattern consisting of 
square cells (SC) or hexagonal cells, respectively. Thus, a natural way to generate a stripe 
filled with the QP pattern in a generic system (e.g., a convective layer), in which a stable 
QP pattern is not available, but stable periodic SC and/or hexagonal patterns do exist, 
is to produce it as a transient layer ("domain wall" H) between two large domains filled 
with the periodic cells, the (half-) sets of the two or three wave vectors in the two domains 
being oriented under the angle, respectively, 45 or 30 degrees relative to each other. The 
pattern in the transient layer will then be a superposition of the two periodic patterns, 
having a full set of the four or six wave vectors, respectively. The so generated stripe will 
feature a QP pattern, provided that its width is essentially larger than the wavelength (a 
size of the elementary cell in the corresponding periodic pattern). It is known that the 
latter condition can be achieved, under special but not unrealistic conditions, for domain 
walls in the nonequilibrium systems @]. 

In fact, it may be simpler to generate, following this way, a twelvefold QP pattern 
between two hexagonal domains, as it is usually much easier to find a stable hexagonal 
structure than a stable SC one j7| . However, the theoretical analysis of the transient layer 
between two SC domains is much simpler, therefore in this work we concentrate on the 
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latter case. In any case, it should be stressed that, although stable square-cell patterns 
are rarer than their counterparts in the form of hexagons or rolls, examples of stable 
square cellular planforms are known in gas flames |J , thermal convection || , and in some 
optical systems [[H||. The interest to the SC planforms has been recently revived by the 
discovery of this pattern in a double-layer Marangoni convection (see the e-print [p] and 
references therein). 

The rest of the paper is organized as follows. In section 2, we formulate the model. 
Some analytical results, which predict, in particular, that the transient layer between the 
two SC patterns has a complex structure, consisting of three layers, two broad ones and 
a narrow sublayer sandwiched between them, are obtained in section 3. In section 4, we 
display results of a direct numerical solution of the stationary real GL equations, which 
comply with the analytical predictions. The paper is concluded by section 5. 

2 The Model 

In this work, we assume a simplest configuration, shown in Fig. [TJ that is going to give 
rise to the QP transient layer. The layer is parallel to the y axis, and the (half) sets of the 
two wave vectors are chosen so that in the left domain the vectors k 12 are parallel to the x 
and y axes, while in the right domain both vectors k3 5 4 have the angle 45 degrees relative 
to the axes. Accordingly, the complex wave field describing the spatial distribution of 
physical variables is assumed to be 

4 

M ( r > t) = J2 B i( r i f ) ex P ( ik i r ) ' ( X ) 

where r is the two-dimensional coordinate, and Bj(r,t) are slowly varying amplitude 
functions. 

In an experiments, the chosen configuration can be created, at least in a part of the 
system, by means of specially selected boundary conditions, which, in turn, are imposed by 
the sidewalls of the experimental cell. Although we do not analyze the sidewall boundary 
conditions in this work, it is obvious that, to support the configuration that we consider, 
one will need to have a large-aspect-ratio cell with two sidewalls forming an angle 45 
degrees, which is quite possible. More general configurations, with different orientations 
of the wave vectors relative to the layer between the two domains, can be considered 
similarly to what is done below, but their technical treatment will be more cumbersome. 

A usual approach to the description of spatially nonuniform patterns of the domain- 
wall type is based on a system of coupled real Ginzburg-Landau (GL) equations for the 
slowly varying amplitudes, assuming that they do not depend on the coordinate y running 
along the stripe (wall), which is a reasonable approximation if the stripe is long enough 
fll"2| , |i~3f . Then, as it was shown in the mentioned works, the effective diffusion coefficient 
(the one in front of the term (Bj) xx ) in each GL equation is 

D ] = {kf)\ (2) 

where ki is the x-component of the vector. For the configuration shown in Fig. [Tl the 
system of the real GL equations can be easily cast into a final form 

A t = (l/2)A xx + A-A 3 -2(g 2 A 2 + g 1 B 2 1 +g 1 Bl)A, (3) 



4 



{B 1 ) t = {B 1 ) xx + B 1 -Bl-2(2g l A 2 + g 2 Bl)B l , (4) 



(B 2 ) t = B 2 -Bl-2(2 gi A 2 + g 2 Bi)B 2 , (5) 

where B\ 2 are the amplitudes corresponding to the wave vectors k 12 in the left domain 
(Fig. p, and, using the obvious symmetry of the configuration shown in Fig. [I], we have 
set £3 = 54 = A. 

Note that the diffusion coefficient in Eq. (H) is zero due to Eq. (fj). Generally speaking, 
in this case one should take into account a higher-order derivative term ~ {B 2 ) xxxx fT2|, |1"3|. 
However, this is not necessary while the diffusion terms do not vanish in the two other 
equations (see details below). 

In Eqs. (^) through (|5|), the linear gain coefficient and the coefficient of the nonlinear 
self- interaction of the spatial mode are normalized to be 1, g\ and g 2 being the coefficients 
of the nonlinear interaction between the modes with the angles a = 45 and 90 degrees be- 
tween their carrier wave vectors. Normally, the nonlinear interaction coefficient decreases 
with the increase of a, so that 

92<9i< 1- (6) 
The necessary and sufficient stability conditions for the SC pattern are well known 

92 < ^ , & > \ (1 + 2g 2 ) , (7) 
while the conditions providing for stability of the eightfold QP planform are 

92 < 2 ,9i < 1 (! + 2^2) • (8) 

An obvious feature of the two sets of the stability conditions is their incompatibility, i.e., 
SC and QP can never be stable simultaneously. 

In fact, the inequality g 2 < 1/2 in the set of the conditions (0) is a cause for the relative 
rarity of stable SC planforms, as, despite the general property (j^), the actual dependence 
g(a) is usually weak, so that g 2 is not essentially smaller than 1. Nevertheless, the full 
set of the SC stability conditions ([!]) can be satisfied in the above-mentioned physical 
systems. 

A consequence of the conditions (^) necessary for the stability of the QP pattern is 
gi < 1/2, which, with regard to Eq. (|]), makes the stability of the QP pattern still less 
feasible than that of the SC one. An objective of this work is to propose a way to produce 
a transient QP pattern between two stable SC domains with different orientations (Fig. 
ID in the case when the uniform QP planform is unstable. Note that a similar approach 
is known as a way to generate of a broad stripe of a SC pattern between two domains 
of orthogonally oriented rolls in the case when the uniform SC pattern is unstable, while 



the rolls are stable 13 



In the case when the SC stability conditions (0) are met, it is quite reasonable to 
assume that the coefficients gi <2 are close to 1/2 (because it is physically implausible to 
have gi^ 2 much smaller than 1/2), which suggests to present them as 

9i,2 = 2^ ~^ 2 )' l^i.2| < 1- (9) 

As it will be seen below, the smallness of /ii^ naturally provides for the strip of the QP 
pattern to be broad, which is exactly the condition justifying the consideration of the 
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transient-layer patterns. Note that // 2 must be positive according to Eq. (0), while \i\ 
may formally have either sign. However, it will be shown below that the necessary solution 
does not exist if \i\ < 0. 

Below, we will also use a parameter 

m = 2//i/// 2 , (10) 

which is, generally, ~ 1. In terms of m, the second stability condition ([?]) takes a very 
simple form, m < 1. 

Eqs. (|3]), (||), and (|5|) must be supplemented by boundary conditions (b.c.) to guar- 
antee that, at x — > ±oo (recall the system is formally assumed to be infinitely large), the 
pattern considered asymptotically coincides with either of the two SC planforms composed 
of the modes Bi and B 2 or B 3 = B 4 = A. Obviously, this implies 

lim A{x) = 1/J2 -fj, 2 = A lim , lim B h2 {x)=0, (11) 
lim Aix) = 0, lim B 12 (x) = A lim . (12) 

X— 00 X— > — 00 ' 

Formally, this set may seem overdetermined, as we add six b.c. to the fourth-order system 
of Eqs. @ and (|4]) (Eq. (||) contains no x-derivatives). However, it is easy to check that 
the seemingly superfluous b.c. for B 2 are nothing else but direct corollaries of the four 
legitimate b.c. for A and B\. 



3 Analytical Results 



A stationary version (d/dt = 0) of Eqs. @, (f|) and (H) can be essentially simplified, as 
in this case Eq. (BT) becomes just an algebraic relation, that has two solutions: B 2 = 0, or 



l-(2- mf x 2 )A 2 -(l-fx 2 )B 2 1: 



(13) 



m being the parameter defined by Eq. (|i0|) . First, we consider the case when the expres- 
sion ( |T3"D holds (obviously, it may hold as long as it yields B\ > 0). Substituting it into 
the stationary versions of Eqs. (|3]) and (|]), we obtain 



A" + mA- (2 (2m - 1) - m> 2 A + (2 — m// 2 ) 



.4 



fx 2 L B'{ + B 1 - (2 - /i 2 ) Bi + (2 - m/i 2 ) A 1 B x 



0. 
0. 



(14) 
(15) 



the prime standing for d/dx. If, instead of Eq. (pT3|), we take _B 2 = 0, the stationary 
equations take the form 



A" + 2A - [2 (2 - /i 2 ) A 2 + (2 - m// 2 ) 5 2 



.4 



+ B x - B\ + (2 - m/x 2 ) A 2 fii 



0. 
0. 



(16) 
(17) 



It is noteworthy that, although the SC patterns may be stable at [ii < 0, i.e., m < 0, 
Eq. (H) cannot have a solution for A(x) exponentially decaying at x — ► —00 (see the b.c. 
(|T2|)) unless m > 0, hence the present problem has no solution with m < 0. 
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Obviously, a solution to Eqs. fll4|) and (|T5|) can satisfy, with regard to Eq. fll3|) , the 



b.c. dl2|) . However, the same set of equations (0) and (|15|) cannot satisfy the b.c. (|TT 
setting i?! = and A = const, one obtains from Eq. (|14"D 



A = A Q = 



2 (2m — 1) — m 2 fi 2 

which is obviously different from the necessary limit value A 2 im = 1/ (2 — fi 2 ). Moreover, 
the value of B 2 corresponding to A\ as per Eq. fllTf ) is different from zero, which also 
violates the b.c. ([□]). 

In fact, the stationary state corresponding to the asymptotic value ([18]) with B 2 7^ 
is another uniform pattern, which is a superposition of three spatial harmonics. This 
pattern is periodic in one spatial direction and quasiperiodic in the other one. As it was 
demonstrated in |Q, this pattern is always dynamically unstable, hence a solution having 
it as an asymptotic state is physically irrelevant. 

However, a solution to Eqs. (|T4|) and (|i~5D makes sense as long as it provides for B 2 > 
according to Eq. (p~3|) . Taking the asymptotic state A 2 = A^, B\ = 0, one finds that it 



gives rise to negative B\ exactly in the case m < 1, which is considered in this work, as 
this is the case when the SC pattern is stable, see above. Thus, Eqs. (|14D and (|lf|) should 
be used in the region —00 < x < x , where, by definition, X — Xq IS cL point at which 
B\ (x), as given by Eq. (|T3|) , vanishes, i.e., 

(2 - m/i 2 ) A 2 (x ) + (1 - /i 2 ) Bl(xo) = 1. (19) 



At the point x = Xq, one must switch from Eqs. ([14]) and ( |T5]) to Eqs. ( |I~6] ) and (|T] 
setting B 2 = at x > Xq. The continuity dictates to take the values A(x — > Xq) and 
Bi(x — > Xq) as the b.c. to Eqs. ( ]T5| ) and ( JH\) at x = Xq, the b.c. at x = +00 being fixed 



by Eqs. ©. 

The transition from Eqs. ( |T4"D and ( |ToD to Eqs. (|TBp and (|T^) at x = Xo provides for 
the continuity of all the functions A(x) and Bi >2 (x), but their first derivatives suffer a 
jump at x = xq. In fact, if the above-mentioned fourth-derivative term is added to Eq. 
(|^), the jump of the derivative will be smoothed down in a narrow boundary layer. It 
should be stressed that the presence of the jump does not violate the applicability of the 
description in terms of the GL equations; the only problem is the absence of a detailed 
description of the narrow boundary layer in which the jump is smoothed by the fourth- 
order derivative term, but details of the inner structure of the narrow layer do not affect 
the global picture. 

Thus, the transient layer between the two SC domains consists of two sublayers, de- 
scribed, respectively, by Eqs. ( |I~4D and (|15|), and by Eqs. (^) and (|T7[) , with the dis- 
continuity of the derivative between them. The first (left) sublayer contains all the four 
spatial harmonics, hence is locally seems as a QP pattern. An important finding is that, 
in the case of small jj, 2 considered here, the width of this sublayer is large, scaling ~ /i 2 , 
according to Eqs. (|T^) and (|T|). The possibility to produce a broad QP stripe justifies all 
the consideration of the transient layer. The second (right) sublayer contains only three 
spatial harmonics, as B 2 = in it, hence it locally looks like a pattern quasiperiodic only 
in one direction H), and its width is ~ 1 (i.e., not specifically large), according to Eqs. 



dH) and (E 

This qualitative analysis of the transient layer's structure will be corroborated and 
illustrated by direct numerical results displayed in the next section. However, before 
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using the numerical methods, one can notice that the structure of the second (right) 
sublayer can be described in an approximate analytical form if we consider a special case 
when the value A(x = xq) at the internal boundary between the two sublayers is already 
close to the asymptotic value A\ im , and, accordingly, the value B\ is small (in other words, 
this is case when < 1-mC 1). Then, we may set 

A(x)=A ]im -a(x), (20) 

where a(x) is positive but small and vanishes at x — +oo. Using the smallness of a(x) 
and Bi, one can simplify Eqs. ( ]TB| ) and (|iT|): 

(1/V2)a" - [V2(2 - /i 2 )a - Bj] - ^B\ = 0, (21) 

B'[ + [V2{2 - fi 2 )a - Bf] B 1 - ~ - 2// 2 ) B 1 = 0. (22) 

The b.c. for a(x) at x = xq can be obtained from the expansion of the exact b.c. ([19]): 

Bf(x ) 

a(lo) = 7i(^)- (23) 



The lowest-order approximate analytical solution to Eqs. ( pT|) and (|22|), satisfying the 
necessary boundary conditions, is very simple: 

Bi(x) = Bi(xq) exp(— Ax), a(x) = a(xo) exp(— 2Xx), (24) 



with A = y (/i2 — 2/ii) /2. Note that the condition necessary for A to be real is again 
exactly tantamount torn < 1. 



4 Numerical Results 

To check the qualitative predictions for the structure of the transient layer obtained in 
the previous section, we performed a two-stage numerical integration of Eqs. ([TJ]) and 
(|T5|), and then of Eqs. (|16|) and ([H])- To this end, the first pair of the equations was 
solved with the b.c. (|T2"D, continuing the solution until it hits a point where it satisfies 
Eq. (|19D (recall the latter condition is equivalent to the vanishing of -B|). Actually, the 
numerical integration was performed, instead of the original coordinate x, in terms of a 
variable £ = tanhx. This transformation is convenient because it maps the semi-infinite 
intervals (±oo,0) of the variable x into finite ones (±1,0). 

The values A(xq) and Bi(xo), obtained from the solution of Eqs. (|i~4"D and ([T5|), where 
then used to find a solution to Eqs. (|l^) and (|T7|) in the interval tanhxo < £ < 1, 



satisfying the b.c. ( |Tl"| ) at £ = 1. In accord with what said above, a condition of the 
continuity of the first derivatives across the point £ = tanhxo was not imposed. 

Three typical examples of the thus obtained numerical solutions are displayed vs. the 
coordinate x in Fig. [| for three characteristic values m = 0.75, m = 0.5, and m = 0.25, 
the parameter /12, that must be small enough, being fixed in all the three cases as \i% — 0.1. 
A characteristic feature clearly seen in all the cases is that, in accord with the prediction 
of the above analysis, a width of the left sublayer is essentially larger than that of the 
right one. It is also noteworthy that the solutions changes very little with a large change 
in m, i.e., the transient layer between the two SC domains is expected to be quite robust. 
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5 Conclusion 



In this work, we have proposed an approach that makes it possible to the generate a 
two-dimensional quasiperiodic pattern in nonlinear dissipative systems where a direct 
generation of stable uniform quasiperiodic planforms is not possible. An eightfold pattern 
can be created in the form of a broad transient stripe between two domains filled by square 
cells, which are oriented under the angle of 45 degrees relative to each other. Using the 
symmetry of the configuration considered, the structure of the pattern was described in 
terms of a system of three coupled real stationary Ginzburg-Landau equations, which were 
analyzed by means of analytical and numerical methods. It was found that the transient 
quasiperiodic pattern exists exactly in a parametric region in which the uniform square- 
cell pattern is stable. Further, it was found that the transient layer consists of two different 
sublayers, with a derivative jump between them (that can be smoothed into an additional 
narrow boundary layer, if higher-order derivatives are added to the Ginzburg-Landau 
equations). The width of the sublayer that features the eightfold quasiperiodic pattern is 
found to be large, while the other sublayer (filled with a less interesting pattern, which is 
quasiperiodic only in one direction) has a width ~ 1. A broad stripe of a twelvefold QP 
pattern can be similarly generated as a transient layer between two domains of hexagonal 
cells oriented at the angle of 30 degrees. 

It still remains to perform simulations of the full time- dependent Ginzburg-Landau 
equations, in order to directly test the dynamical stability of the broad transient layers. 
However, the numerically found robustness of the layers against the variation of the crucial 
control parameter m suggests that they have a good chance to be dynamically stable. 
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Figure Captions 



Fig. 1. The configuration giving rise to the transient layer filled with the eightfold quasi- 
periodic pattern between two domains of square cells oriented under the angle 45 degrees 
relative to each other.. 

Fig. 2. Numerically found structure of the transient layer corresponding to the con- 
figuration shown in Fig. 1. The small parameter /12 defined by Eq. (|9]) is fixed to be 0.1, 
while the control parameter m, defined by Eq. (|i~0|), takes values 0.75 (a), 0.50 (b), and 
0.25 (c) (recall only the values < m < 1 make sense in the present context). 
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